Magnetoelectric bistabilities in ferromagnetic resonant tunneling structures 
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The conditions for the occurrence of pronounced magnetoelectric bistabilities in tlie resonant 
tunneling through a ferromagnetic quantum well are theoretically investigated. The bistability 
appears due to the mutual feedback of the carriers Coulomb interaction and the carriers exchange 
coupling with magnetic impurities in the well. It is shown that the well Curie temperature depends 
strongly on the relative alignment of the quantum well level and the reservoirs chemical potentials, 
which can be modified electrically. Switching between a "current-on/magnetism-off" and a "current- 
off/magnetism-on" mode becomes possible, if the well temperature lies in-between the bistable values 
of the well Curie temperature. 
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In ultimate magnetoelectric devices the magnetic prop- 
erties should be ideally controllable to a vast extent by 
external bias or gate fields. For this purpose, band- 
engineered magnetic resonant tunneling structures are 
very promising, since they exhibit a rich variety of tun- 
able magneto-transport properties [l|. Especially, the 
impetuous development of novel dilute magnetic semi- 
conductors (DMSs) 0, y, 3 in the last decades, which 
are made magnetic by randomly doping with transition 
metal elements, e.g., by incorporating Mn in a GaAs 
crystal host, has considerably enriched the possibilities 
of growing different magnetic semiconductor heterostruc- 
ture systems. In DMSs the ferromagnetism can depend 
strongly on the actual particle density, which has been 
confirmed in several experiments, in which ferromag- 
netism has been generated by tailoring the particle den- 
sity by electrical or optical means [5|, |6| . 

In magnetic resonant tunneling structures made of 
para- or ferromagnetic DMSs even small energetic spin 
splittings of the well subbands can become observable 
in the transport characteristics 0, d, 3 • Based on their 
spin-dependent transmission magnetic resonant tunnel- 
ing structures have been proposed for realizing efficient 
spin valves and spin filtering devices [l|, |lO| , or for dig- 
ital magnetoresistance jlll Il2l |. The magnetic proper- 
ties of ferromagnetic quantum wells made of DMSs are 
well described in the framework of a mean field model 
[13l . UA . [iSl . [iq . which reveals that the Curie tempera- 
ture of the well depends on (i) the 2D-spin susceptibility 
of the carriers and (ii) on the overlap of the subband 
wave function with the magnetic impurity density pro- 
file. Both parameters should be in principle tuneable by 
the applied bias, which would provide a purely electri- 
cal control of the ferromagnetism in magnetic quantum 
wells. 

In conventional nonmagnetic resonant tunneling diodes 
(RTDs) it is well known, that an intrinsic hysteresis in 
the negative-differential-resistance (NDR) region of the 
current- voltage (IV) characteristics can occur [13]. This 
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FIG. 1: (Color online) (a) Schematic scheme of the band pro- 
file of the magnetic double barrier structure. The exchange 
interaction of the magnetic ions is mediated by the carriers 
tunneling in and out of the well, (b) Equivalent circuit model 
of the resonant tunneling structure introducing the emitter 
and collector capacitances Ce,Cc and resistances Re, Re, re- 
spectively. 



bistability of the tunneling current has been explained 
to result from the nonlinear feedback of Coulomb inter- 
action of the stored well charge [l^, [l^. In magnetic 
RTDs this naturally suggests the possibility of hysteretic 
magnetic states, as has been predicted in J20ll2l|. 

In this article a detailed study of possible magneto- 
electric bistabilities in magnetic RTDs is provided. The 
carriers dynamics is described by a self-consistent sequen- 
tial tunneling model, which includes the feedback effects 
of both the carriers Coulomb interaction and the mag- 
netic exchange coupling with the magnetic ions. The 
model yields a simple expression for the steady state 2D- 
spin susceptibility, which allows to calculate the critical 
temperature Tc of the quantum well depending on the 
applied bias and the relative alignment of the quantum 
well level with respect to the chemical potentials of the 
emitter and collector reservoirs. If the well is operated 
at a temperature, which lies between the bistable values 
of the well Curie temperature, the magnetic RTD can 
be switched between a "current-on/magnetism-off" and 
a "currcnt-off/magnctism-on" mode. 



The band profile of a generic double-barrier resonant 
tunneling structure with a ferromagnetic quantum well 
made of a DMS, e.g., of GaMnAs, is sketched in Fig.[lja). 
The vertical transport through the structure can be de- 
scribed by a sequential tunneling model, since the high 
density of magnetic impurities in the well will likely cause 
decoherence processes. By using the transfer Hamilto- 
nian formalism a Pauli master equation for the statisti- 
cal distribution of the particles in the well can be derived 
Uli [23] ■ In the case that only a single resonant level E^ 
resides in the energy window of interest, which is defined 
by the difference of the emitter's and collector's chemi- 
cal potentials, simple rate equations for the spin-resolved 
well particle densities Na{t) with (cr =t, i) are obtained 



diV, 
dt 



T,iE,) Ne,a + T,{E^) TV,,, - r{E,) N„. (1) 



Here, N^^^^,.^ are the densities of particles with the res- 
onant longitudinal energy E^ in the emitter (e) and col- 
lector (c) reservoir, respectively. The energy-dependent 
tunneling rates rre,,-i,r = Fe -I- F, can be calculated 
by Bardeen's formula [23|, which essentially evaluates 
the overlap of the lead and well wave functions in the 
barriers. For high barriers these tunneling rates be- 
come proportional to the longitudinal momentum pz of 
the incident particles [23], i.e., Fe^c oc {E^Y^"^ with E^ 
denoting the longitudinal energy. By assuming that 
the particle reservoirs are described by Fermi-Dirac dis- 
tributions the particle densities are given by Ni_a = 
DQkBT\n{\+eyiYi[{^ii- Ea)/kBT]}, i = (e,c), with 
Dq = mj^'KTi} is the two-dimensional density of states 
per spin for carriers with the effective mass m, ks de- 
notes Boltzmanns' constant, T is the lead temperature, 
and [ii are the emitter and collector chemical potentials 
with {ic = lie — eT4ppi where 14ppi is the applied bias. 

In the framework of a mean field model an analytic 
expression for the steady state exchange splitting A of 
the well level can be derived [U, [ij, [l^, [iq 
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where Jp^ is the coupling strength between the impurity 
spin and the carrier spin density (in case of GaMnAs p- 
like holes couple to the d-like impurity electrons), z is the 



longitudinal (growth) direction of the structure. 



.W 



is the impurity density profile, i\jq(z) labels the well wave 
function, and s = 1/2 is the particles spin. The Bril- 
louin function of order S is denoted by Bg^ where S is 
the impurity spin, which for Mn equals 5/2. By consider- 
ing a homogenous impurity distribution A is effectively 
determined by the voltage dependent spin polarization 

? = s(^T-^i)- 

The nonlinear feedback of the Coulomb interaction of 

the well charges is approximately taken into account by 

calculating the electrostatic well potential in terms of 
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FIG. 2: (Color online) (a) Contour plot of the well Curie 
temperature T^ (K) as a function of the well level position 
_B,„ and the applied bias, (b) Schematic illustration of the 
different occupation probabilities of the quantum well level 
from the reservoirs for the regions A-D, as indicated in the 
contour plot (a). The emitter and collector band edges are 
denoted by _Eoe and i5oc, respectively. 



an equivalent circuit model of the resonant tunneling 
diode, as shown in Fig. [^b), where the capacitances Ce 
and Cc are determined by the geometrical dimensions 
of the barriers and the well [221 . The potential results in 
U:^ = [e2(iV-7Vback)-Cceyappi]/C, where iV = TVt + iV^, 
C = Ce-I-Cc, e denotes the elementary charge, and iVback 
is the positive background charge in the well, which origi- 
nates from the magnetic donors. Since the actual position 
of the quantum well levels E^ = Eq + 11^ — crA depends 
on both the magnetic exchange splitting A and the elec- 
trostatic potential C/„, all equations become nonlinearly 
coupled, making a selfconsistent numerical solution nec- 
essary. 

In order to find criterions for the occurrence of mag- 
netic bistabilities and to interpret the numerical results 
in the following it is very useful to study the dependence 
of the well Curie temperature T, on the applied bias and 
the well level position. The mean field model yields an 
analytic expression for the collective Curie temperature 
of a magnetic quantum well 



knTr = 
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where the two-dimensional spin susceptibility is defined 
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FIG. 3: (Color online) The (a) quantum well level position 
-Eu,, (b) current j, (c) Curie temperature Tc, and (d) well 
splitting A as a function of the applied bias. The solid lines 
indicate the voltage up-sweep values (F), whereas the dashed 
lines correspond to the voltage down-sweep values (R) . In (a) 
the i5„-voltage curves are embedded in the contour plot of the 
well Curie temperature Tc (K) and in (d) the solid red line 
corresponds to the actual well temperature Tweii = 4.2 K. 
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Within the introduced sequential tunneling model 
Eq. ([T]) the steady state spin susceptibility simplifies to 
X2DiE) = ~s{dNo/dE), where iVo = (re7Ve+r,A^c)/r is 
the steady state solution of the rate equations H]). Hence, 
the dimensionless susceptibility can be written as 



X2d{E) 
sDa 



2— e,c 



FD 



N^_d_ 

DodE 



(5) 



with fpD,i = (e,c) denoting the Fermi-Dirac function 
for the emitter and collector reservoir, respectively. This 
allows to calculate the Curie temperature Tc as a func- 
tion of the applied voltage (note that Tc depends via /p^ 
explicitly on the voltage, since fi^ ~ ^J■e — eVappi) and 
the well level position, as displayed in Fig.[2)^a). For the 
simulations I used generic parameters corresponding to 
a GaMnAs well: m = 0.5 mo, £r = 12.9, d^ ^ dc = 20 
A, w = 10 A, /ie = 70 meV, nimp = 1.5 x lO^Ocm"^^ 
^pd — 0.06 cV nm'^, where de,dc and w are the emit- 
ter barrier, collector barrier and quantum well widths, 
TTiQ denotes the free electron mass, and e^ is the rel- 
ative permittivity of the well. The background charge 
JT-back = 0.1 njmp is Considered to be only of about 10% 
of the nominal Mn doping density [2J| and the lattice 
temperature is set to T = 4.2 K. 

The Curie temperature contour plot can be divided 
into four qualitatively different regions A-D. which are 



characterized by different probabilities for occupying the 
quantum well level from the reservoirs, as schematically 
illustrated in Fig.[2Ub). In region A, for instance, the well 
level can be occupied by particles originating from both 
reservoirs. The Curie temperatures Tc ' in regions A 
and C differ roughly by a factor 2 compared to T^^ of 
region B. This sudden change of Tc can be explained 
as follows: by assuming energy-independent tunneling 
rates and nearly symmetric barriers, i.e., Fg k, Vc the 
dimensionless spin susceptibility of Eq. (O simplifies for 
region A to x^ = 1/2(/|q + /^p) w 1, whereas in the 
other regions one obtains: x^ = /fd/2 « 1/2,%'^ = 
/fd ^ ^^ ''^^'^ X^ ~ 0- This simple estimation, hence, 
yields the desired result TJ^-'^V^f « 2. 

These differences in the Curie-temperatures of the var- 
ious regions can now be exploited to realize hysteretic 
magnetoelectric states. According to the nonlinear feed- 
back of the stored well charge, the resonant level i?^, and 
the IV-characteristic show a hysteretic behavior, as dis- 
played in Fig. [3]^a) and (b). For the up-sweep (F) of 
the applied bias the well is charged before the E^^ be- 
comes off-resonant, i.e., when it drops below the emit- 
ter band edge, whereas for the voltage down-sweep (R) 
the well is almost uncharged before E^j becomes reso- 
nant again. This leads to different self-consistent elec- 
trostatic potentials for up- and down-sweeping voltages, 
explaining the occurrence of the intrinsic bistability. If 
the hysteresis of E^, now switches exactly between the 
Tc-regions B and C, as it is the case in Fig. El^a), then 
also the voltage-dependent Curie temperature will ex- 
hibit a pronounced hysteresis, as shown in Fig.[3{c). The 
electric hysteresis will then be accompanied by a mag- 
netic hysteresis if the actual lattice temperature T of the 
quantum well, which is displayed as straight solid line 
in Fig. ^c), fulfills the condition TJ^ < T < Tf . This 
is illustrated in Fig. [3ljd) : as long as the resonant level 
stays in region B the well is nonmagnetic (A = 0, since 
T > Tc) but when £"„, enters the region C the well be- 
comes immediately magnetic (A ^ 0). Also notice, that 
at low voltages the well is always magnetic. At roughly 
Kippi ~ 30 mV the well becomes nonmagnetic, since 
Ew crosses the boundary between the regions A and B, 
which provides a purely electrical control of the well mag- 
netism. As a whole, the magnetic well switches between 
a " current-on/magnetism-off " state for the up-sweep and 
a "current-off/magnetism-on "-state for the down-sweep 
of the applied voltage. Moreover, the Tc contour plot 
in Fig. mja) also suggest the possibility for realizing the 
switching between a "currcnt-on/magnctism-on" and a 
" current-off/magnetism-off " mode in the case that the 
hysteresis of E^, switches between region B and D and if 
0<T<Tf. 

In summary, I have shown by using a selfconsistent se- 
quential tunneling model that the Curie temperature of 
a magnetic quantum wells strongly depends on the rela- 
tive alignment of the well level and the reservoirs chemi- 
cal potentials, which can be modified by external bias or 
gate fields. Magnetoelectric bistabilities become possible 



if the hysteresis of the weh level position E^u switches 
between regions of different well Curie temperatures and 
if the actual well temperature lies in-between these two 
values. 
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